UNCLASSIFIED 

414749 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


ROTICK:  When  governnent  or  other  drawings,  sped* 
fleet  Ion#  or  other  dete  ere  used  for  eny  purpose 
other  than  In  connection  with  a  definitely  related 
governaent  procuraewnt  operation,  the  U.  8. 
Oovernaent  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever)  and  the  fact  that  the  Oovern¬ 
aent  eay  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  laplication  or  other¬ 
wise  as  in  any  aanner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri£its 
or  permission  to  aanufhcture,  use  or  sell  any 
patented  invention  that  nay  in  any  way  be  related 
thereto. 


UNIVERSITY  OF 


MARYLAND 


0> 

z> 

H 


THE  INSTITUTE  FOR  FLUID  DYNAMICS 

and 


APPLIED  MATHEMATICS  jul  r 


Technical  Note  HI- 32 3  June  1963 


SUftMANY  OF  "THEORJE  ANALYTICJUE  DES  PB0BLEKE8 
STOCHASTIQUES  KELAXIFS  A  UN  GROUFE  DE  IlOREB 
TELEFHOJflQUES  AVEC  DXSFOSITXF  D'AHTHTIE"  • 


Felix  Pollaczek 
University  of  Maryland 
College  Part,  Maryland 


•  Supported  in  part  by  the  Office  of  Naval  Research-MOMR  39$  (17)  sad 
the  U.  S.  Any  Research  Office  (Durhaa)  Grant  DA-ARO  (D) -31-12^-0  13A 


Introduction 

The  purpose  of  this  memoir  is  to  treat  in  the  utmost  general¬ 
ity  the  most  obvious  queueing  problems  concerning  a  bundle  of  s  fully 
available  telephone  trunks  (or  a  system  of  s  equally  accessible  counters) 
where  calls  (or  customers)  are  treated  according  to  the  rule  "first  come, 
first  served". 

Although  being  concerned  with  probability  problems,  we  don't 
borrow  from  probability  calculus  anything  but  the  very  notions  of  proba¬ 
bility,  of  dependent  r-  independent  chance  variables  and  of  distribution 
functions  (d.  f.)  of  one  or  several  variables.  As  our  methods  are  exclu¬ 
sively  analytical  ones  .us  employ  widely  the  elementary  parts  of  the  theory 
of  analytical  functions  of  one  complex  variable. 

In  general  we  make  the  following  assumptions; 

For  all  holding  times  Tn  and  interarrival  intervals  (of  success¬ 
ive  calls)  Yn  (n*0,l,2,... )  the  probabilities 

Psidr  (  Tm  <  i)  (t)  Frur‘x  (  Y*.  < 


d 


are  arbitrarily  given  distribution  functions,  the  only  restrictive 
hypothesis  being  the  condition 

E(t)  -j*  }  oiW  <  -o. 


As  initial  condition  we  suppose  that  from  the  instant  of 
generation  of  a  given  cal*  ( which  is  denoted  as  call  0"  )  on. 

If 

there  elapse  respectively  jtol  |  0  ,  *  *  •  £  O 

units  of  time,  before  all  s  trunks  are  free  from  previously  generat¬ 
ed  calls. 


Then  one  of  our  basic  results  runs  as  follows.  Let  X.\ 
be  the  waiting  time  of  the  n  th  call,  q  °) 

a  complex  parameter  and  <y  (i*i  <0  a  complex 

variable.  She  generating  function  of  the  mathematical  expectations 

Ee‘ft~  i*  «• 

5  <>.■>)  "S.  * 

is  given  by  the  following  expression: 

vm*S 


(1)  }<n>  -YM  *-&j£  + 

~im*f  'i++t  *  * 

♦-  +£?  1”  r  e  V"** 


-im+S  •£•>+£ 


(f  »o) 
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where  the  functions  VJ  are  given  by  the  system  of  s  linear  simultane¬ 
ous  integral  equations 


(2) 


e,(t)-i  i 
* 


i-  £*  C~  ^  V  “  0  ^X+I 


(\  =  0jt)  •••,  A-/) 


and  the  equation 


Z  S  i>) 

X*«  i*.  , 


*  O  . 


Here 


has  to  be  extended  over  all 


combinations  of  the 


s  indices  /, 
fined  by 


A. 


£i(y)  €i(y)  are  respectively  de- 


£*(>)  *1  df.(i) 

O 


(<->,*  Rty*0) 


while 


l-o 


The  functions 


are  analytic  for 


Vi(  Wi'h*  >) 

R(^.)  >  o 

R  ( <V»)  >  o 
R  ( ^)  >o 


and.  synnetric  in  <y, 

(l)  these  functions  appear 


In  all  integral  formulae  similar  to 


only  in  the  form 


The  system  of  integral  equations  (2),  (3)  can  he  solved  by 
a  finite  number  of  steps  in  the  following  two  cases: 

l.°  For  £,(y)  rational,  i.  e.,  when  the  d.  f.  4}("k)  is  of  the 

form 

f,(t)  «  i  -  £ 

where  Py( i)  is  a  polynomial  in  t,  fi(^)  (and  therefore  ) 


being  an  arbitrary  function 


a 


j 


2.o  For  £,(y)  of  the  form  i.  e.  when  ft(t)  is 
lattice  function  with  distance  being  rational. 

As  is  known  by  the  theorem  of  Paul  Lfevy,  which  in  the  present 
case,  owing  to  the  fact  that  T»  io } can  be  written  in  a  simpler  form, 
the  dLf  fjt)  ~  Rt^  fa* can  derived  from  E  C~*‘**' 


by  the  formula 
(*) 


Lormuxa  ^ ^ 

.  _ i  f  r  (t>o&>0). 

f~«>  771  \  e  Ee  r 


-i*/+s 


More  generally  generating  functions  like 

t,  >fl~  EVfWw 


eao 


•5-  « 

2-L>i  1  Y  c 


#  «•  A 
^  *0 


are  always  given  by  formulae  of  the  form  (1),  (2),  (3);  only  the  ritfxt 
sides  of  the  integral  equations  (2)  axe  in  every  particular  case  to  be 
modified  according  to  the  problem  treated. 

If  i„*0;  }  f  o  t  that  is,  if  at  the  instant  of  gmsra- 

tion  of  call  "0"  all  trunks  axe  unoccupied,  all  terms  on  the  right  side 
of  (l),  save  V/o) ,  disappear,  ao  that  than  this  formula  talma 


u 


the  simpler  form  J  ^  V0(°). 

While  all  reasoning  on  this  memoir  is  conceived  in  such  a  way  as  is  re¬ 
quired  for  treating  the  many  server  problem  .  I  ,  our 

formulae  are  valid  also  for  s  =  1  . 

The  methods  employed  here  apply  also  to  all  probability  prob¬ 
lems  arising  for  a  fully  available  bundle  of  trunks  without  waiting  de¬ 
vice,  i.  e . ,  to  the  calculation  of  various  loss-of -calls  probabilities. 

CHAPTER  I.  Construction  of  the  Laplace-Stieltjes  transforms 
of  different  distribution  functions  of  waiting  times. 

In  order  to  calculate  the  mathematical  expectation  E  C 
which,  by  the  help  of  (4)  ,  yields  the  d.f.  fmli)  we  proceed  as  fol¬ 
low*.  Let  Xfr  (/*■*  be  the  instants  of  generation  of  our  calls, 

so  that  Xm-h  *“  X*v  =  Ym.  ('***>■>"■)  “d 

Xm,  +  f  ,  Xa.  +  -Xns. 

be  the  last  ■  ends  of  conversations  asked  for  before  time 

(and  consequently  having  indices  <  m)>  taken  in  an  arbitrary  order. 

As  there  are  s  trunks,  the  waiting  time  t*.  of  the  n  th  call  will  be 
/m**v  (,4m  ***,4.a  )  ,  if  this  number  is  >  O  ,  and  zero  in  the 

opposite  case.  Writing 

CL+  — 
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we  have  therefore 
(5)  %L 


4 

t rru/r u 


We  now  use  the  formula 

(«,  ^rw'v'-A)] = '-s^H 

((?(,!  >o)  C'  4 

where  the  integration  paths  C />  *  * ' ,  £*.  are  parallels  to  the  imgi- 

nary  axis,  situated  to  the  right  of  this  axis,  and  traversed  from  be¬ 
low  to  above. 


This  formula  shall  be  proven  by  induction.  Supposing  its 
validity  for  s  -  1,  we  shall  show  that  it  is  also  valid  for  s  ;  as  it 
is  valid  for  s-1,  it  is  true  for  all  s. 

For  s  ■  1,  (5)  runs  as  follows: 

e*'**  -  '  -  ITT  f  e‘*  T^r  -f  ( *>•,«&  • 


The  integration  path  shall  be  shifted,  while  rwsining  to  it¬ 

self,  indefinitely  towards  the  right  for  4^0,  and  towards  the  left 
for  a  >0  .  Bus,  we  see  that  for  <14  0  this  integral  is  aero,  but 
for  a>0  it  is  eqpal  to  the  sub  of  residues  in  t^m0  and  f}m  “f  , 
that  is  to  I  *  C  ^  ,  and  this  proves  (5)  for  s  ■  1. 

For  A  >  I  ,  we  can  adait  that  CL^  *  <t,9  3*,  *%»-/  • 


8 


As  multifold  integrals  of  this  type  are  absolutely  convergent,  which 
is  shown  in  supplement  I,  we  can  treat  (5)  as  an  iterated  integral. 
Proceeding  in  the  same  manner  as  previously,  we  find  that 

+  t) 

(*>>0)  ■ 


Thus  for  0  formula  (5)  is  evidently  right.  For 
the  contribution  of  the  second  tern  to  the  right  to  the  remaining 
(.  -  1)  low  integral  U  0  because  In  11  eapooMt. 

<3y-  are  by  hypothesis  C  O  ,  while  in  virtue  of  the  first 
term 


_  .  —  the  remaining  (a  -  0  fold  integral  has  the 
r  t*  *9* 

form,  with^S-/)  instead  of  s,  as  the  integral  (5) •  Thus,  because 
for  ~  0+)  > 

we  have  transformed  (5)  into  the  analogous  formula  for  / 
which  is  true  by  hypothesis.  Therefore,  (5)  is  proven  for  all  ^ • 
In  virtue  of  (5)»  we  get 

e  (*c  l  t+rV 
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This  formula  shows  In  an  evident  manner  that 
tion  of  the  . 


e~u 


a  symmetric  func- 


We  want  to  calculate 


E  =  E(eTf^\A oj. 


The  waiting  time  5*.  is  a  definite  function  of  the  chance  variables 
%  ,  ~Q-i  (holding  times)  and  Yi,  'C-i  (interarrival  intervals 

between  successive  calls )  and  of  the  given  parameters  -ie, }  '*'  >  ^oa 
which  describe  the  initial  conditions  at  the  Instant  of  generation  X« 
of  call  0".  therefore,  replacing  in  all  integral  formulae  Ty  and 

II 

lay  the  corresponding  snail  letters,  we  have 

^  mO  g| 

(t)  E(e**l*J  “j  "  e'**" 

when  e  But  be  expressed  tv  the  Jiwy  ec=erdi«  to  (6.) 

It  la  easily  seen  that  the  quantities  are  equal,  save 
for  their  order,  to  the  quantities 


<8)  ‘*v  * l- ' u  ~r->'">zz 


m*l 


<  <  /mi*.**  £  "•  -S  • 


where  /mm*v  5  /r»ww. 


As  the  maters  ,  which  are  relative  to  the  calls  of 

indices  <  /*»“/ ,  depend  only  on  the  chance  variables  ^ 

but  not  an  7i*-i  and  ,  we  have 
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-at*.  , 

vtiex*  C  is  given  by  (6).  V/e  substitute  for  the  the  expres¬ 

sions^),  which  shows  that  (9)  is  a  synmetric  function  of  the 
and  can  then  perform  the  integrations  indicated  in  (9)  under  the  in¬ 


tegral  signs  by  replacing  factors  of  the  form  £ 
respectively  by 


or  e 


il-,* 


do) 


After  these  operations  we  represent  (9)  as  a  sun  of  Fourier 
integrals  in  the  real  variables  t  _>  which  is  done  by  virtue  of 

different  trans format ions  explained  in  supplement  III.  In  the  sees 
ganger  we  can  calculate  the  conditional  natbesntical  expectation 

m  m  v 


and  represent  it  as  a  ra  of  Fourier  integrals  in  the  variables 


11 


Repeating  this  process  fa- &  times,  we  get  finally  for  Ewntj 
which  is  the  9tt.e It jes -Laplace  transform,  under  the  given  initial 
conditions  jto9}  of  fa(t)  *=*  (t*.  <  t)  an  expression  of  the  fol¬ 

lowing  form 


V  -J—tO 

(•*0  lit  4 1 

+ i  j(|Ze>V>Su.W^ 


(U) 


•i-*6 


+  ••  • 


-WM* 


<«W 

where  j  slgilfles  ihat  the  path  of  integration  is  parallel  to, 
->♦# 

end  situated,  at  the  right  of,  the  1  magi  nary  axis  of  the  complex  plane. 


12 


The  functions  VMX  are  defined  by  the  initial  conditions  and 
the  recurrence  formulae 

(>.•0,1,  •*  >  +~0  ) 
JMO 


* 


(12) 


Pm  lug 


.JK  J'h’V)  =  +  jfe  f-5^1  • 

/|A  -!•* 

(\~o, j 

ir.  'Z  Z  'jktyrfai*  )• 

*  *■*  ,;  ..(V«i 

Vx(>.,  ".>»i>)  =  £  b’V, 


/»v 


multiplying  ill)  by  ^  ,  and  the  last  two  formulas  respectively 

by  /y,#'  and  and  sunning,  we  obtain  the  formulae  contained  in 
the  introduction.  These  formulae  reduce  the  calculation  of  fm(*) 
essentially  to  the  resolution  of  the  system  of  equations  (2),  (3), 
and  this  problem  is  treated  in  chapters  IV  and  V.  In  a  manner  anal¬ 
ogous  to  (k),  the  composite  probabilities 

(*,  *'  U.)  -  ^  ft*  *«  w  **' 1  *») 

can  be  expressed  by  meant  of  the  mathematical  expectations 
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Replacing  in  (11)  ^  ,  hi  and  /'*■  respectively  by  J  )  4,*  and  ^ 
we  have 

=  */*> + W  K4,  •*•*)  t’ 


(13) 


-'(**0 


.!«>+• 


In  accordance  with  (13),  we  oust  now  multiply  this  expres¬ 
sion  by  e  (eq.  (6a))  which  is  also  a  Fburier  integral  in  the 

same  variables  and  represent  this  product  as  a  sum  of  Fourier 

integrals  in  the  -4«  v  •  This  can  be  done  by  means  of  formula  (3.25) 
of  suppl.  III.  Next  we  have  to  carry  out  the  operation 
E  (  Ji0y)(see(7)  )  which  yields  on  expression  of  the  form  (11) 
with  recurrence  formulae  of  the  fora  (12),  the  only  difference  being 
that  the  initially  given  lTc^  are  replaced  by  other  expressions. 

Likewise  the  generating  function 

2  yy*'  i*») 

/*«  O,  m.'mo 


is  given  by  a  sum  of  integrals  of  the  form  ( 1),  with  integrands 
satisfying  ( 3 )  and  a  system  of  integral  equations  of  the  form  (2), 
although  with  different  right  sides. 


By  the  same  method  more  general  composite  probabilities 
such  as  Pji^  (%.<*,  <i"  \ 

can  also  be  reduced  to  the  resolution  of  a  system  of  integral 
equations  of  the  form  (2),  (3). 

Chapter  II.  Construction  of  the  generating  functions 

of  different  probabilities . 

Vika. 

Denoting  byAthe  probability  that,  at  the  instant  of  generation 
of  the  call,  exactly  a  calls  ( generated  after  tine  X*  ) 

be  waiting,  the  generating  functio.  ^  .  It*  a 

AM  9 

is  represented  by  means  of  the  function  *<* •?) 

defined  and  calculated  previously. 

The  generating  function  Z-  Zj  *  'h  «.  of  the 
probabilitie^that,  at  the  instant  of  production  of  the  n  tta  call, 
exactly  X  =  &  trunks  be  occupied  is  also  calculated.  This 

problem  can  be  reduced  to  the  resolution  of  a  system  of  linear  equ¬ 


ations  differing  from  (2)  only  by  their  right  sides  and  with  y--d$ 

i"  ”** 

replaced  by  JUm.  J  di , 
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CHAPTER  III .  Problems  concerning  the  distribution  of  the 
Markovian  Parameters .llieory  of  the  Phenomenon  of  Temporary  Blocking. 

The  quantities  defined  in  the  beginning  of  Chapter  I 

can  be  considered  as  'Markovian  parameters  "  relative  to  the  instant 
X/k)  because  they  characterize  exhaustively  all  that  which  the 
knowledge  of  events  concerning  calls  generated  before  the  instant  X*. 
can  tell  us  about  phenomena  which  are  posterior  to  X*. . 

We  have  seen  previously  that  min*  Is  the  waiting  time 
of  the  n  bn  call,  and  have  explained  our  method  how  to  obtain 
the  distribution  function  off  “5*.  More  generally  we 

(*)*f  i 

can  study  the  distributions  of  the  quantities 
which  are  the  symbolic  "waiting  times1'  elapsing  between  X*  and 
the  Instant  when  at  least  /U  trunks  of  the  bundle  become  free  of 
conversations  asked  for  before  time  Xr*. . 


By  means  of  our  method  we  reduce  the  calculation  of  the  generating 
functions  of  the  mathematical  expectations  (m.  e.) 

ECeT*  I  0* 9 

and, 

E(e~*  ( 1(p)>°,#fJ>0) 


to  the  problem  of  resolution  of  a  system  of  integral  equations  anal¬ 
ogous  to  (2),  (3). 


l6 


In  order  to  obtain  the  joint  distribution  of  all  s  quanti¬ 
ties  /mi-*,  (p  * '>">*)  >  it  is  necessary  to  calculate  the  m.  e. 


We  do  this  by  means  of  a  class  of  three-indices -operators  T£r 
which  we  have  used  already  to  treat  other  problems  of  this  kind,  and 
for  the  particularities  of  which  we  refer  to  a  previous  paper  (Applica¬ 
tion  d'  operateurs  integro-  combinatoires  dans  la  thdorie  des  lntdgrales 
multiples  de  Dirichlet",  Ann.  Inst.  H.  Poincard,  v.  11,  19^9, 
p.  113-133).  For  the  calculation  of  the  generating  function  of  the 
last  mathematical  expectation  the  formulae  (1),  (2),  (3)  remain  still 
validj however,  the  right  sides  of  (2)  must  be  replaced  by  certain 
functions  of  •  We  treat  then  the  particular  case 

i(t)  «  l-c *  j  =  IT^T  > 

in  which  these  functions  have  a  particularly  simple  font. 

In  the  second  part  of  this  chapter  we  generalise  the  assump¬ 
tions  which  were  made  up  to  now,  in  order  to  be  able  to  take  into  con¬ 
sideration  the  so  called  phenomenon  of  temporary  blocking.  This  pheno¬ 
menon  due  to  technical  reasons  has  the  following  effect.  From  the  in¬ 
stant  when  the  n  th  call  is  assigned  to  a  non -occupied  trunk,  this 
one  and  all  other  non -occupied  trunks  are  blocked  during  a  certain  time  I*. 


IT 


We  consider  the  quantities  U/n.  as  stochastic  variables,  stochastical¬ 
ly  independent  of  all  other  variables  save  ~C  and  assume 

that 

Hob  (x.  <  t,  ^<0  )  =  , 

ilt.e )  denoting  a  giver,  distribution  function  of  two  variables 

which  is  independent  oi  n  .  is  xn  the  preceding  chapters  we  suppose 


that  oc 

E(t)  -  J  i  <  <*> . 

0 


By  T^,  we.  de6i©mte  now  the  waiting  time  imposed  on  the  n  th  call 
toy  the  sole  effect  of  previous  calls  (of  indices  <  /n,  ).  To  T m. 
must  toe  added  the  blocking  delay  ,  so  that  the  n  th  convex-  * 

sat  ion  begins  at  time  ©m. . 

The  calculation  of  the  £T  ^  ^  I  )  la  saaewhat  aore 
conplicated  than  under  our  previous  assumptions.  For  the  generating 
function  of  these  mathematical  expectations  we  obtain  again  an  expres¬ 
sion  of  the  for*  (1),  tout  the  functions  Vx  ere  now  determined  I ay  a 
system  of  lipear  equations,  in  the  left  sides  of  which  there  appear 
simple  aa  well  as  double  definite  integrals,  while  the  rlgrt  sides  are 
«P*1  to  i+y,  . 


0 


Under  the-  particular  assumptions  0  ^  ^ 

which  implies  that  fi*  and  -.re.  stochastically  independent,  and 

f,(i>  =  ,  ■'<-  ha"1,  emoioyed  the  above  mentioned  system 

of  integral  cuuatiour  to  '•ale  liato  t  h<*  quantity  V.fo) 

,  to  which 

Z  V  E(e~,T~ \i.t) 

( in  Reduction 

probability  d'  attente  au 


reduces  itself  t:>r  j£,(  £  c  ^ 
de  differ  nts  piutidi:*:;-  ov;c 


±..  S  O 

f 


tdldphine,  d  ia  resolution  i«.  . vstdoes  't«  Equations  integrales". 


Ann.  Inst.  H.  io incur'7. 
CHAPTKR  IV. 


ft-,’.  P-  ijt  -  l'j  (p.  165). 


-t 


Pecclutlor,  of  the  iru  ”rc.i  equat ior.s  (2),  (3)  for  = 
lhat  is .  for  £t($)  -  ~y  . 

When  trie  Let  Mac  -Stioitjes  transform 
is  a  rational  fmctMii,  ttn  s  integral  equations  (2)  can  be  trans¬ 
formed  into  functional  equations,  for  a ny  fjt)  and  its  tranform 

Position  Li,’  ...ho  study  o:  r  he  rose  of  an  arbitrary  rational 
t.hb  -  ha  ve  treated  ir.  this  er.apv.r  the  case  of  the  s  implest 
rational  E,(f)  i.  <'  of  €/$)  —  j _  j  in  a  more  thorough 
manner.  Notice  here,  thrt  nr.:  ring  •*  |  — ^  ,  ue  have 

CM  *{uf,(t)  “  »  which  signifies  that  we  have 

taken  the  mean  holding  tins  as  unity  of  times. 
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Now  put  €,  ( t)  —  ,  _  j  in  the  integral  equations  (2).  By 
virtue  of  the  properties  of  the  analytic  functions 
explained  in  supplement  IV  the  integrands  appearing  in  (2)  have  then 
in  the  right  half  plane  of  the  complex  variable  ^  a  unique  pole 
Jwj  and  are  OOtl”*)  for  R(O>0  ,  It  I  h  00 
Accordingly  the  integrals  in  (2)  are  equal  to  the  respective  residues 
of  their  integrands  at  f  *  /  ,  so  that  these  equations  are  trans¬ 


formed  into  t 

d4)  -  =  f  +  r 

(x*  A- 0. 


As  has  already  been  mentioned,  in  the  simplest  case  studied  in  queue¬ 
ing  theory,  that  is,  when  at  the  instant  X,  of  generation  of  the  ini¬ 
tial  call  0"  all  s  trunks  are  non-occupied  so  that  $0."’ 

ft 

equation  (1)  becoams 


the rs fore,  we  are  principally  Interested  in  the  construction  of  V#fy) 
and  V.(0). 
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To  that  end  we  have  to  put  in  the  equation  (14)  i  ■*  '  =  ^ 

(that  is,  to  take  all  r^9  equal  to  the  unique  pole  of  &t(t)  ) 

For  the  s  +  1  functions 

(15)  =  (x*0, 

we  obtain  in  this  manner  from  (14)  and  (3)  the  s  +  1  linear  non 
homogeneous  equations 

D -»"-»>]  W  -^r  [(*r>) - 

(16)  -  «,(-!-<)  V„,  (**•)]  =  (X'°>''-‘r,), 

X** 

Therefore,  V.(^)  can  be  expressed  as  the  quotient  of  tw>  func¬ 
tions  of  y  ,  numerator  depending  Unwarily  on  the  s  still 
unknown  quantities  V>)  and  the  denominator  containing  the 
factor 

(17)  sa-y,  . 

Posing  «  X  in  (16)  we  get  for  the  the  s  -  1 

homogeneous  linear  equations 

(18)  [i-^cj-xVIV*)  -Tr^-x-OVk/atO  «  o  (w  ,r'\ 

It  is  easily  seen,  as  a  consequence  of  touch* 's  theorem, 
that  for  1*^1  <  *  the  function  (17)  has  in  the  right  half  plane  of 
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the  variable  y  a  unique  zero  ^  (iy)  .  Since  all  functions 
^')***j  y*  i  *y)  *  6,1(1  hence  the  functions  (15),  are  finite  for 

,  •  •  •  ,  Rt'yO  1  o  )  R l^)  10  ,  the  aforementioned  numerator 

of  V0t^)  must  also  disappear  for  =  and  this  property  yields 

an  s  th  linear  relation  which  together  with  equations  (18)  permits 
us  to  calculate  all  s  quantities  V\  ( M  .  For  V0  ( 6)  ve 
obtain  thus  the  formula 

a,)  =  “VI 

n\~C  * 


where 


is  given  by  the  equation 


X*o  * 

>  v^>]. 
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From  (19)  we  obtain  for  the  generating  function  of  the  distribu¬ 
tion  functions  ~  ft*,  (i  I  °)  by  virtue  of  equation  (4) 


(21) 


Fty'f  r ~ 

ff\-0 


——  f  1  + 

'-If  L 


One  can  obtain  thus  the  distribution  functions  ^  ( i) 
by  developing  the  right  side  in  a  Taylor  series,  and  in  accordance  with 
the  present  assumption  Jtc)  $  0^  •*  *  ,  ^  0  one  sees  that  in  parti¬ 
cular  y*(i)  =1  . 

In  order  to  find  the  limit  distribution  function 

<22>  f<*>  =  ?J*) 

we  represent  ^ ( t)  as  a  complex  integral,  i.  e., 

f-1  *»  =  iii  Ff,>’ 

where  K  denotes  a  little  circle  the  center  of  which  is  ^  *  0  • 

If  f*C«  is  such  that 

£»  =  /, 

there  exists  a  constant  >  |  such  that  for  l^|  <  rf  the  zero 

fljofy)  of  the  expression  (17)  is  a  holomorphic  function;  the  sane  is 
therefore  true  for  P«(^)  (equ.  (20)  and  (21)  ). 
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Therefore,  if  has  no  zero  on  the  circumference 

1^1  -  I  1  which  is  true  if  the  limit  (22)  exists,  we  find,  extending 
K  beyond  the  unit  circle,  that 

_i 

fU)  =  l+fl.  0)e 

and  can  moreover  establish  an  asymptotic  evaluation  for  the  difference 
fit)  -  ||VN  It)  .  Nov  the  existence  (with  a  single  exception)  of 

the  limit  (22)  has  been  proven,  for  arbitrary  -fi(t)  and  TA ( t) 
by  J.  Kiefer  and  J.  Wolfowitz.  But  here  we  need  not  resort  to  this 
theorem,  for  we  prove  that  in  the  present  case 

so  that  the  existence  of  the  limit  (22)  is  evident. 

Also  for  fc’td-fjt)  $  ~  } 

whence  f(f)  —  0  ,  it  is  possible  to  establish  asymptotic  formulae 

for  fj±)  . 

The  end  of  this  chapter  contains  the  construction^  for 
£i(  $)m  |'_  j  }  of  the  function  V,  ( ^ ,  j ^ 0  which  appears  in  (1), 
as  well  as  several  applications  of  the  last  formulae  to  the  problems 
of  Chapter  II. 
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CHAPTER  V.  Resolution  of  the  integral  equations  (2),  (3) 
under  different  assumptions  concerning  the  Laplace -st ielt Jes  transforms 
£  1  C  0  and  S)  * 

The  first  of  the  two  cases,  mentioned  in  the  introduction, 
in  which  the  equations  (2),  (3)  can  be  resolved  by  means  of  a  finite 
number  of  steps,  is  outlined  under  the  assumption 

f,(i)  =£  *,+**+-+** 

i-l 

that  is  for 

(23)  £,(s)  =  Z,  * 


Proceeding  in  the  same  manner  as  on  the  occasion  of  the  de¬ 
duction  of  equation  (l4),  we  can  replace  then  all  integrals  in  (2)  by 
the  sums  of  the  respective  residues  of  their  integrands  in  $  * 

Thus  these  integral  equations  are  transformed  into  functional  equations, 
the  method  of  solution  of  which  is  explained  in  detail  for  s  »  2  and 
outlined  for  arbitrary  s  .  In  order  to  construct  Vt(^ 
for  s  -  2  ,  it  is  necessary  to  calculate  a  certain  determinant 
^iven  here  for  n  =  2  and  n  =  3  )  and  to  determine  those  ') 

roots  =  "•  of  the  eqjjation 
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which  for  sufficiently  small  !^|  are  situated  in  the  right  jy-half- 
plane. 


When  the  numbers  ^y0)  are  all  different,  the  limit  distri¬ 
bution  function  ii/v*  is  of  the  form 


i*  -»<*> 


f 


It) 


=  l  -3^e',sW±  (t>a, **J, '*$■*■). 


S>  =  | 


For  every  s  ,  the  V*  If)  appearing  in  (2),  (3)  are,  under 
the  present  assumptions,  rational  functions  of  ^  a®*1  . 

The  second  case  in  which  the  equations  (2),  (3)  can  be  resol¬ 
ved  by  means  of  a  finite  number  of  operations,  is  studied  under  the  as¬ 
sumptions 


The  method  of  resolution  of  (2),  (3)  is  explained  in  detail  for  •  »  2 
and  then  outlined  for  arbitrary  s  ;  in  every  case  the  are  ration¬ 

al  functions  of  y.  In  order  to  construct  V0(<^) 
fbr  s  =  2  ,  it  is  necessary  to  calculate  those  <m(fe+0  zeros 
of  another  determinant  function  (given  here  for  n  «  1,  2,  3) 


which  for  sufficiently  small  Ir^l  are  situated  in  the  left  y  -  halfplane. 

Finally  we  suppose  that  from  the  initial  instant  X0  an  enormous 
number  of  calls  is  generated.  Then,  of  course,  JLs**-  Pm  (t)  ~  0  (**«) 

m-#oo  J 

but  we  nevertheless  study  this  case  in  order  to  be  able  to  establish  for 
asymptotic  expressions  for  m  — *  «>»  . 

Cur  suppositions  ear  best  be  rendered  by  assuming  that 

(25)  fz(i)  =  0  -I  ( t  (yy  Jthd  -  I. 

Then  the  form  of  the  integral  equations  ( 2)  as  well  as  the  theory  of 


their  resolution  under  both  Hypotheses  for  C t  (f^)  treated  in  this  chap¬ 
ter  is  considerably  simplified,  ie  rind  that,  for  an  £,(f)  according 
to  (23)  as  well  as  under  the  hypotheses  (?t) ,  our  integral  equations 
can  be  resolved  by  resolving  several  systems  of  linear  algebraic  equa¬ 
tions,  without  having  to  resort,  to  certain  solutions  of  transcendental 
equations . 


In  order  not  to  interrupt  too  often  the  exposition  of  our 
theory,  several  demonstrations  have  been  placed  at  the  end  of  this  memoir 
(supplements  SI  -  S7). 


